The theory of magnetic field induced domain-wall propagation in magnetic nanowires 
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A global picture of magnetic domain wall (DW) propagation in a nanowire driven by a magnetic 
field is obtained: A static DW cannot exist in a homogeneous magnetic nanowire when an external 
magnetic field is applied. Thus, a DW must vary with time under a static magnetic field. A moving 
DW must dissipate energy due to the Gilbert damping. As a result, the wire has to release its 
Zeeman energy through the DW propagation along the field direction. The DW propagation speed 
is proportional to the energy dissipation rate that is determined by the DW structure. An oscillatory 
DW motion, either the precession around the wire axis or the breath of DW width, should lead to 
the speed oscillation. 



Magnetic domain-wall (DW) propagation in a 
nanowire due to a magnetic field 0, S, H, II H reveals 
many interesting behaviors of magnetization dynamics. 
For a tail-to-tail (TT) DW or a head-to-head (HH) DW 
(shown in Fig. 1) in a nanowire with its easy-axis along 
the wire axis, the DW will propagate in the wire un- 
der an external magnetic field parallel to the wire axis. 
The propagation speed v of the DW depends on the field 
strength^, Q. There exists a so-called Walker's break- 
down field Hyy @ • v is proportional to the external field 
H for H < Hw and H 3> Hw- The linear regimes 
are characterized by the DW mobility \i = v/H. Ex- 
periments showed that v is sensitive to both DW struc- 
tures and wire widthP,i,Q. DW velocity v decreases as 
the field increases between the two linear H-dependent 
regimes, leading to the so-called negative differential mo- 
bility phenomenon. For H ^S> Hw, the DW velocity, 
whose time-average is linear in H, oscillates in fact with 
time 0,0]. 
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FIG. 1: Schematic diagram of a HH DW of width A in a 
magnetic nanowire of cross-section A. The wire consists of 
three phases, two domains and one DW. The magnetization 
in domains I and II is along +z-direction (9 — 0) and -z- 
direction (6 — n), respectively. Ill is the DW region whose 
magnetization structure could be very complicate. H is an 
external field along +z-direction. 

It has been known for more than fifty years that 
the magnetization dynamics is govern by the Landau- 
Lifshitz- Gilbert (LLG)Q equation that is nonlinear 
and can only be solved analytically for some special 
problems @, [§]. The field induced domain- wall (DW) 



propagation in a strictly one-dimensional wire has also 
been known for more than thirty years [fj, but its exper- 
imental realization in nanowires was only achieved[l|, Q, 
& S, HI m recent years when we are capable of fabrica ting 
various nano structures. Although much progress 0, llfj| 
has been made in understanding field-induced DW mo- 
tion, it is still a formidable task to evaluate the DW 
propagation speed in a realistic magnetic nanowire even 
when the DW structure is obtained from various means 
like OOMMF simulator and/or other numerical software 
packages. A global picture about why and how a DW 
propagates in a magnetic nanowire is still lacking. 

In this report, we present a theory that reveals the 
origin of DW propagation. Firstly, we shall show that 
no static HH (TT) DW is allowed in a homogeneous 
nanowire in the presence of an external magnetic field. 
Secondly, energy conservation requires that the dissi- 
pated energy must come from the energy decrease of the 
wire. Thus, the origin of DW propagation is as follows. 
A HH (TT) DW must move under an external field along 
the wire. The moving DW must dissipate energy because 
of various damping mechanisms. The energy loss should 
be supplied by the Zeeman energy released from the DW 
propagation. This consideration leads to a general re- 
lationship between DW propagation speed and the DW 
structure. It is clear that DW speed is proportional to 
the energy dissipation rate, and one needs to find a way 
to enhance the energy dissipation in order to increase the 
propagation speed. Furthermore, the present theory at- 
tributes a DW velocity oscillation for H ^> Hw to the 
periodic motion of the DW, either the precession of the 
DW or oscillation of the DW width. 

In a magnetic material, magnetic domains are formed 
in order to minimize the stray field energy. A DW that 
separates two domains is defined by the balance between 
the exchange energy and the magnetic anisotropy energy. 
The stray field plays little role in a DW structure. To 
describe a HH DW in a magnetic nanowire, let us con- 
sider a wire with its easy-axis along the wire axis (the 
shape anisotropy dominates other magnetic anisotropics 
and makes the easy-axis along the wire when the wire is 
small enough) which is chosen as the z-axis as illustrated 
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in Fig. 1. Since the magnitude of the magnetization M 
does not change in the LLG equation^, the magnetic 
state of the wire can be conveniently described by the 
polar angle 8(x, t) (angle between M and the z-axis) and 
the azimuthal angle <p(x,t). The magnetization energy 
is mainly from the exchange energy and the magnetic 
anisotropy because the stray field energy is negligible in 
this case. The wire energy can be written in general as 



E = J F(6,(t),V6,\7<j))d 3 x, 

F = f(6,<t>) + ^[(V6») 2 + sin 2 6»(V0) 2 ] - MH cos 0, 



(1) 



where / is the energy density due to all kinds of magnetic 
anisotropics which has two equal minima at 8 — and 
t = Oj </ ) ) = = 7r >'^)); J~ term is the exchange 
energy, M is the magnitude of magnetization, and H is 
the external magnetic field along z-axis. In the absence 
of H, a HH static DW that separates 8 = domain and 
8 = tt domain (Fig. 1) can exist in the wire. 

Non-existence of a static HH (TT) DW in a magnetic 
field-In order to show that no intrinsic static HH DW 
is allowed in the presence of an external field (H ^ 0), 
one only needs to show that following equations have no 
solution with 6 = at far left and 8 = tt at far right, 
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= JV 2 HM sin 8 - J sin 8 cos 8 '(V '4>f = 0, 
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= 0. 



(2) 



Multiply the first equation by V# and the second equa- 
tion by V(/>, then add up the two equations. One can 
show a tensor T satisfying V • T = with 

T =[/ - HMcosO + ^(|V6»| 2 + sin 2 8\ V0| 2 )]1- 

J(V(9V(9 + sin 2 8V(tN<j)), 

where 1 is 3 x 3 unit matrix. A dyadic product (V8V8 
and V0V0) between the gradient vectors is assumed in 
T. If a HH DW exists with 8 = in the far left and 
8 = tt in the far right, then it requires — /(0, <j>) + HM = 
-/(tt, <p)-HM that holds only for H = since /(0, (f>) = 
f(n,(j)). In other words, a DW in a nanowire under an 
external field must be time dependent that could be ei- 
ther a local motion or a propagation along the wire. It 
should be clear that the above argument is only true for 
a HH DW in a homogeneous wire, but not valid with de- 
fect pinning that changes Eq. |2]). Static DWs exist in 
fact in the presence of a weak field in reality because of 
pinning. 

What is the consequence of the non-existence of a 
static DW? Generally speaking, a physical system un- 
der a constant driving force will first try a fixed point 
solution ll| . It goes to other types of more complicated 



solutions if a fixed point solution is not possible. It means 
that a DW has to move when an external magnetic field 
is applied to the DW along the nanowire as shown in 
Fig. 1. It is well known \10\ that a moving magnetiza- 
tion must dissipate its energy to its environments with a 
rate, = sM. f+™ (dm/dt) 2 d 3 x, where m is the unit 

vector of M, a and 7 are the Gilbert damping constant 
and gyromagnetic ratio, respectively. Following the simi- 
lar method in Reference 12 for a Stoner particle, one can 
also show that the energy dissipation rate of a DW is 
related to the DW structure as 



dE 

~dt ~ ~ (1 + a 2 )M 



M x H, 



eff 



d 3 i 



(3) 



where H e f f = 



SF 
SM 



is the effective field. In regions I and 

II or inside a static DW, M is parallel to H e ff. Thus no 
energy dissipation is possible there. The energy dissipa- 
tion can only occur in the DW region when M is not 
parallel to H e f / . 

DW propagation and energy dissipation-For a magnetic 
nanowire in a static magnetic field, the dissipated energy 
must come from the magnetic energy released from the 
DW propagation. The total energy of the wire equals 
the sum of the energies of regions I, II, and III (Fig. 1), 
E = Ei + En + Em. Ei increases while En decreases 
when the DW propagates from left to the right along the 
wire. The net energy change of region I plus II due to 
the DW propagation is 



d(Ei + En) 
dt 



-2HMvA, 



(4) 



where v is the DW propagating speed, and A is the cross 
section of the wire. This is the released Zeeman energy 
stored in the wire. The energy of region III should not 
change much because the DW width A is defined by the 
balance of exchange energy and magnetic anisotropy, and 
is usually order of 10 ~ lOOnm. A DW cannot absorb 
or release too much energy, and can at most adjust tem- 
porarily energy dissipation rate. In other words, dE d ' t " 
is either zero or fluctuates between positive and nega- 
tive values with zero time-average. Since energy release 
from the magnetic wire should be equal to the energy 
dissipated (to the environment), one has 
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cry 



dt 



(1 + a 2 )M Jin 



M x H e ff ) d A x. 

(5) 



or 



2(1 + a^)HA Jin 



(rh x H e f A d 3 x+- 
\ J 1 ) ; 



1 dE 



in 



2HMA dt ' 
(6) 

Velocity oscillation-Eq. ([6]) is our central result that 
relates the DW velocity to the DW structure. Obviously, 
the right side of this equation is fully determined by the 
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DW structure. A DW can have two possible types of 
motion under an external magnetic field. One is that a 
DW behaves like a rigid body propagating along the wire. 
This case occurs often at small field, and it is the basic as- 
sumption in Slonczewski modelQ and Walker's solution 
for H < Hw- Obviously, both energy-dissipation and 
DW energy is time- independent, dE d ' t " = 0. Thus, and 
the DW velocity should be a constant. The other case 
is that the DW structure varies with time. For example, 
the DW may precess around the wire axis and/or the DW 
width may breathe periodically. One should expect both 
dE d ' t " and energy dissipation rate oscillate with time. Ac- 
cording to Eq. ([6]), DW velocity will also oscillate. DW 
velocity should oscillate periodically if only one type of 
DW motion (precession or DW breathing) presents, but 
it could be very irregular if both motions are present and 
the ratio of their periods is irrational. Indeed, this os- 
cillation was observed in a recent experiment [3( . How 
can one understand the wire-width dependence of the 
DW velocity? According to Eq. ©, the velocity is a 
functional of DW structure which is very sensitive to the 
wire width. For a very narrow wire, only transverse DW 
is possible while a vortex DW is preferred for a wide wire 
(large than DW width). Different vortexes yield differ- 
ent values of |m x H e f / 1 , which in turn results in different 
DW propagation speed. 
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FIG. 2: The time-averaged DW propagation speed versus 
the applied magnetic field for a biaxial magnetic nanowire 
of cross section 4nm x 20nm. The wire parameters are 
A"i = K 2 = 10 5 J/m 3 , J = 4. x 10" 11 J/m, M = 10 6 A/m, 
and a = 0.1. Cross are for the calculated velocities from Eq. 
0, and the open circles are for the simulated average veloci- 
ties. The dashed straight line is the fit to the small H < Hw 
results, and solid curve is the fit to a(H — Ho) 2 /H + b/H. 
Insets: the instantaneous DW speed calculated from Eq. © 
for H = 50Oe < H w (left) and H = lOOOOe > H w (right). 



Time averaged velocity is 



a 7 



(rh x H e ff) d 3 x, (7) 



2(1 + a 2 )HA J in 



where bar denotes time average. It says that the averaged 
velocity is proportional to the energy dissipation rate. In 
order to show that both Eqs. (O and ([7]) are useful in 
evaluating the DW propagation speed from a DW struc- 



ture. We use OOMMF package to find the DW struc- 
tures and then use Eq. J7J) to obtain the average velocity. 
Figure 2 is the comparison of such calculated velocities 
(cross) and numerical simulation (open circles with their 
error bars smaller than the symbol sizes) for a magnetic 
nanowire of cross-section dimension 4nm x 20nm with 
a biaxial magnetic anisotropy / = — ^-M 2 + ^-M 2 . 
The system parameters are K\ — K 2 — 10 5 J/m 3 , 
J = 4. x 1CT 11 J/m, M = 10 6 A/m, and a = 0.1. The 
good overlap between the cross and open circles confirm 
the correctness of Eq. (J7|) . The v — H curve for H > Hw 
can be fit well by aA(H — H ) 2 /H + b/H (see discus- 
sion later). The insets are instantaneous DW propaga- 
tion velocities for both H < Hw and H > Hw, by Eq. 
(J6j) from the instantaneous DW structures obtained from 
OOMMF. The left inset is the instantaneous DW speed 
at H = 50O e < Hw, reaching its steady value in about 
Ins. The right inset is the instantaneous DW speed at 
H = lOOOOe > Hw, showing clearly an oscillation. They 
confirm that the theory is capable of capturing all the 
features of DW propagation. 

The right side of Eq. ([7]) is positive and non-zero 
since a time dependent DW requires rh x H e ff ^ 0, 
implying a zero intrinsic critical field for DW propa- 
gation. If the DW keep its static structure, then the 
first term in the right side of Eq. ([6|) shall be pro- 
portional to a A AH 2 , where a is a numerical number 
of order of 1 that depends on material parameters and 
the DW structure. This is because the effective field 
due to the exchange energy and magnetic anisotropy 
is parallel to M, and does not contribute to the en- 
ergy dissipation. Thus, in this case, v 



'1+a 3 



H 



with /u = CL jj^p . Consider the Walker's ID modelf 



in which / = -^M 2 cos 2 9 + ^-M 2 sin 2 6 cos 2 (f>, here 
K\ and K 2 describe the easy and hard axes, respec- 
tively. From Walker's trial function of a DW of width 



A, In tan 



1 

A(t) 



z — J v(r)dT and (f>(z,t) — 4>(t), 



one has (from Eq. (J3])) the energy dissipation rate 



dE 
~dt 



2c*7,4A 
1 + a 2 



[K 2 M 3 sin 2 cj) cos 2 + H 2 M] , (8) 



and DW energy change rate is 

Substituting Eqs. (jHJ) and ^ into Eq. ©, one can easily 
reproduce Walker's DW velocity expression for both H < 
Hw and 3> Hw- For example, for H < Hw — aK 2 M/2 
and A = const., Eq. ^ gives 



aA7 
1 + a 2 



K 2 M sin</>cos< 
H 



H. 



(10) 



This velocity expression is the same as that of the 
Slonczewski model|9l for a onc-dimcnsional wire. In 
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Walker's analysis, 4> IS fixed by K2 and H through 
KiM sine/) cos <j> = Using this <p in the above ve- 
locity expression, Walker's mobility coefficient fx = 2 — 
is recovered. This inverse damping relation is from the 
particular potential landscape in (^-direction. One should 
expect different result if the shape of the potential land- 
scape is changed. Thus, this expression should not be 
used to extract the damping constant [l|, 

A DW may precess around the wire axis as well as 
be substantially distorted from its static structure when 
H > Hw as it was revealed in Walker's analysis. Ac- 
cording to the minimum energy dissipation principle [l3|. 
a DW will arrange itself as much as possible to satisfy 
Eq. (2). Thus, the distortion is expected to absorb part 
of H. The precession motion shall induce an effective 
field g{4>) in the transverse direction, where g depends 
on the magnetic anisotropy in the transverse direction. 
One may expect rh x H e ff — {H — Hq) sin 9z + sin 9 g((j>)y, 
where Ho is the DW distortion absorbed part of H . Using 
\mxH e ff\ 2 = (H-H ) 2 sin 2 9 + g 2 sin 2 9in'Eq. Q, the 
DW propagating speed takes the following h-dependence, 
v = aajA(H-Ho) 2 / H(l+a 2 )+bajA/[H(l+a% linear 
in both A and H for H 3> Hq, but a smaller DW mobil- 
ity. This field-dependence is supported by the excellent 
fit in Fig. 2 for H > Hw- The reasoning a gree s also 
with the minimum energy dissipation principle [13j since 
m x H e f f \ — H sin 9 when M for H — is used, and any 
modification of M should only make |m x H e ff\ smaller. 
The smaller mobility at H 3> Hw , Hq leads naturally to 
a negative differential mobility between H < Hw and 
H ^> Hw^- In other words, the negative differential mo- 
bility is due to the transition of the DW from a high 
energy dissipation structure to a lower one. This picture 
tells us that one should try to make a DW capable of 
dissipating as much energy as possible if one wants to 
achieve a high DW velocity. This is very different from 
what people would believe from Walker's special mobility 
formula of inverse proportion of the damping constant. 
To increase the energy dissipation, one may try to re- 
duce defects and surface roughness. The reason is, by 
minimum energy dissipation principle, that defects are 
extra freedoms to lower \rh x H e ff \ because, in the worst 
case, defects will not change \rh x H e ff\ when M without 
defects are used. 

The correctness of our central result Eq. ((6|) depends 
only on the LLG equation, the general energy expres- 
sion of Eq. |T]), and the fact that a static magnetic field 
can be neither an energy source nor an energy sink of a 
system. It does not depend on the details of a DW struc- 
ture as long as the DW propagation is induced by a static 
magnetic field. In this sense, our result is very general 
and robust, and it is applicable to an arbitrary magnetic 
wire. However, it cannot be applied to a time-dependent 
field or the current-induced DW propagation, at least not 



directly. Also, it may be interesting to emphasize that 
there is no inertial in the DW motion within LLG de- 
scription since this equation contains only the first order 
time derivative. Thus, there is no concept of mass in this 
formulation. 

In conclusion, a global view of the field-induced DW 
propagation is provided, and the importance of energy 
dissipation in the DW propagation is revealed. A gen- 
eral relationship between the DW velocity and the DW 
structure is obtained. The result says: no damping, no 
DW propagation along a magnetic wire. It is shown that 
the intrinsic critical field for a HH DW is zero. This 
zero intrinsic critical field is related to the absence of a 
static HH or a TT DW in a magnetic field parallel to 
the nanowire. Thus, a non-zero critical field can only 
come from the pinning of defects or surface roughness. 
The observed negative differential mobility is due to the 
transition of a DW from a high energy dissipation struc- 
ture to a low energy dissipation structure. Furthermore, 
the DW velocity oscillation is attributed to either the 
DW precession around wire axis or from the DW width 
oscillation. 
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